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MATERIAL  NONLINEARITY 

I'.  Kupptisamy ,  A.  Narnia ,  and  .) .  N.  Krddv 
ABSTARCT 

A  fully  three-dimensional  analysis  of  laminated  plates  with  mater¬ 
ial  nonlinearity  is  presented.  The  modified  Rombe r-Osgood  relation  is 
used  to  compute  the  principal  elastic  moduli  in  the  plane  of  the  laminae. 
The  elastic-plastic  model  with  Hill's  criterion  is  also  used  in  the  study. 
A  shear  deformation  plate  theory  is  also  used  to  compare  the  results 
obtained  by  the  three-dimensional  model. 


INTRODUCTION 

Considerable  research  effort  has  gone  in  the  past  into  the  analysis 
of  laminated  composite  plates.  Much  of  the  earlier  research  was  confined 
to  analyses  based  on  the  laminate  plate  theory  of  Re  issuer  and  Stavsky  [l]. 
Extension  of  the  classical  laminate  theory  to  thick  plates  is  due  to  Yang, 
Norris,  and  Stavsky  [2].  Closed-form  solutions  of  this  theory  were  pre¬ 
sented  by  Whitney  and  Pagano[3],  and  Reddy  and  Chao  [a].  Reddy  and  his 
colleagues  [5-8]  and  Spilker[9]  presented  finite-element  analyses  based 
on  the  shear  deformation  plate  theory.  A  state  of  plane  stress  was  assu¬ 
med  in  all  plate  theories,  and  the  individual  lamina  and  hence  the  laminate 
is  assumed  to  behave  in  linear  elastic  fashion.  The  finite-element  analysis 
of  laminated  plates  using  the  three-dimensional  elasticity  theory  has  been 
conducted  by  Lin  [lO],  Dana  [ll],  Dana  and  Barker  [ 1 2 ] ,  and  Kuppusamy  and 
Reddy  [l3].  The  studies  in  [10-12]  were  limited  to  linear  elastic  analyses 
while  that  in  [ 1 3 ]  was  geometricall  nonlinear.  None  of  these  studies  consi¬ 
dered  material  nonlinearity  or  elasto-plast ic  behavior  in  the  analyses. 
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('riff in,  Kamat  and  Herakovich  [1M  included  material  nonlinearity 
through  a  elastic-plastic  model  with  strain  hardening  to  study 
laminated  plates  subjected  to  inplane  loads.  Pifko.  Levine,  and  Armen 
[15]  presented  a  simiLar  study.  Nonlinear  material  properties  were 
introduced  via  one-dimensional  Ramberg-Osgood  representations  by  Renieri 
and  Herakovich  [16].  All  of  these  studies  were  mainly  intended  to 
compare  laboratory  test  data  in  one-dimensional  tension  tests  with 

those  predicted  hv  FKM/Ramberg-Osgood  relation.  Plate  bending  problems 
were  not  studied. 

The  present  study  involves  a  fully  three-dimensional  analysis  of 
composite  plates  using  the  finite  element  method  and  includes  material 
nonlinearity.  Apparently,  the  present  study  is  the  first  one  to 
consider  the  material  nonlinearity  in  three-dimensional  finite-element 
analysis  of  composite  plates.  The  material  nonlinearity  is  introduced  by, 
1)  uncoupled  one-dimensional  modified  Ramberg-Osgood  model;  and  2) 
elastic-plastic  model  using  modified  Hill’s  criterion  for  anisotropic 
media.  The  finite  elements  developed  in  [5]  and  [13]  are  extended  to 
account  for  material  nonl inearit ies .  Some  of  the  results  from  the 
three-dimensional  nonlinear  analysis  are  compared  with  those  of  shear 
deformable  ( 2— D)  theory. 


FINITE  ELEMENT  FORMULATION 


Here  we  present  a  displacement  finite  element  model  for  an  elastic 
body.  The  kinematic  description  of  an  elastic  body  yields  the  following 
equations  of  equilibrium  (in  the  absence  of  body  forces  and 
moments) 


where 


o..  -  0  (i,j  =  1,2,3) 

*  J 

are  the  stress  components,  and  are 


(1) 

the  Cartesian  coordinates. 


The  linear  stress-strain  constitutive  relationship  i r 
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where  C...  .  contains  the  material  constants,  s  .  is  the  strain  tensot  . 
ljkl  kl 

The  field  equations  are  to  be  adjoined  by  the  boundarv  conditions 
specified  for  the  problem.  The  finite  element  model  associated  with  the 
variational  form  of  the  equation  (see  [17])  over  an  arbitrarv  element 
of  the  finite  element  mesh  is  given  by 


r  f 

[6u.  .  a..]  dx.  dx0  dx~  =  /  t  u  ds 
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where  6  denotes  the  variational  symbol  and  t  is  the  surface  traction  on 

J  m 


The  displacements  u^  are  assumed  in  terms  of  nodal  displacements 
through  interpolation  functions.  These  are  expressed  in  matrix  form  as: 


{u  }  -  [  N 1  (q 


where  q  ■  is  the  nodal  displacement  vector  and  ( N ]  is  the  interpolation 


For  a  three-dimensional  eight-node  brick  element,  one  has 


Nt  -  i  (i  +  ssi)  (i  +  yy±)  (l  +  tt±)  i  «  l * ... 8  (5) 


where  s,  y,and  t  are  the  local  coordinates  (i.e.,  in  the  element). 


The  strain-displacement  relations  can  be  expressed  in  terms  of  the 
column  of  nodal  displacements  as 


[k]  =  l  B I  |q  l 


4 

where  [  B  ]  is  the  matrix  consisting  of  the  derivatives  of  the 
interpolation  functions,  N^.  Substituting  Fqs.  (A  )  and  (6)  into  Fq .  ( ) 
and  taking  the  variation,  one  obtains  the  element  equation 

[Kel  {q  )  =  {R°}  (7) 

where 


[Kel  = 


;b]  [cl  [ B 1 dv 


(8") 


and 


{R  }  =  applied  force  vector=  J  {t  } [N ]  ds 

r  e 


(9) 


The  element  equations  are  assembled  in  the  usual  manner  to  obtain 
the  global  equations,  which  have  the  general  form 


[ K 1  {Q}  =  {R>  (10) 

Equation  (10)  is  solved  for  (Q 1  after  introducing  the  essential 
boundary  conditions.  The  secondary  quantities,  i.e.,  strains  and 
stresses  in  each  element  are  then  computed. 


MATERIAL  MODELS 


Two  different  nonlinear  material  models  are  considered.  They  are, 
1)  Modified  Ramberg-Osgood  model;  and  2)  elastic-plastic  model  with 
modified  Hill’s  criterion  for  anisotropic  media. 

Modified  Ramberg-Osgood  Model 

For  an  anisotropic  linear  elastic  medium,  the  stress-strain 


relationship  (2)  can  be  written  in  a  matrix  form  as: 
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where  A^  to  A^  are  given  in  terms  of  engineering  constants  F^  ,  F^ ,  F^, 

^12’  J21’  ^13’  V31’  "23’  J,i2'  S'*’  (’13  and  G23  ^see  If  the 

stress-strain  relationship  is  nonlinear,  then  the  coefficients  A^_.  vary 

during  deformation.  If  it  is  assumed  that  the  nonlinearity  can  be 
represented  by  uncoupled  stress-strain  behavior  for  each  of  the  strain 
components,  then  a  suitable  curve  fitted  to  the  experimentally  observed 
stress-strain  behavior  can  serve  as  the  material  constitutive  law.  The 
Ramberg-Osgood  relation  used  here  is  one  such  procedure. 

For  each  of  the  strain  components  the  stress-strain  relationship 
can  be  simulated  by  the  Roraberg-Osgood  relation  (see  [19,20])  as 


♦  • 
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where  €  =  strain,  o  -  stress,  E.  =  the  initial  modulus, and  A  and  m  are  the 
parameters  defining  the  curve  (see  Figure  1). 

A  modified  stress-strain  relationship  can  be  derived  from  the  above 
law  as: 

E.  -  E 

0  = - -5  -  -V  IT-iT,-  *  E_  (13) 
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The  tangent  modulus  can  now  be  defined  as 


E  =  ^ 
t  d  f 


K .  -  K 
.  1 _ JL_- 

E.  -  E  m  l+m 

i.  i  1  Pi  V  m 


(14) 
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The  above  equation  is  used  to  calculate  the  F  values  at  every 
strain  level  and  update  the  stillness  matrix  in  the  finite  element 
f ormu 1  at  ion . 

Elastic-Plastic  Model 

An  incremental  formulation  is  used  here.  It  is  assumed  that 


‘7e  °p 


:  ;  i  =  1,2,. ..,6 


(151 


i )  Oe 

where  •  -  incremental  total  strain,  e.  =  incremental  elastic  strain  and 

i  i 

??  -  incremental  plastic  strain. 

i 

The  Hooke’s  (2)  law  in  incremental  form  is  given  by 


o 

=  A 


i,j  =  1,2,. ..,6 


(16) 


where  -  the  constitutive  matrix  given  in  Fq.(ll). 

As  per  the  modified  Hill’s  criterion  for  an  orthotropic  medium,  the 
following  yield  function  is  used  [26]: 


2f(cr.)  =  F  (o2 


+  G(<J3  -  oi)2  +  H(oi  -  a2)2 

+  2M (  °^ ) 2  +  2N(°4)2  =  1 


+ 


2L(  °j)  2 


(17) 


where  t(  •'  )  =  yield  function,  and  material  parameters  F,  G,  H,  L,  M,  and  N 
are  defined  in  terms  of  yield  strengths  as: 


(18) 


A 
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where  X,V,  and  7.  are  the  normal  yield  strengths  in  the  material  principal 
directions  and  R,S  and  T  are  the  shear  yield  strengths. 

The  normality  condition  is: 


i 


Ho  . 
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i  =  1,2. - ,6 


(19) 


where  X  can  be  evaluated  from  the  consistency  condition  give  by 


dl  -  o.  =  0  i  -  1,2 . 6  (20) 

.  1 

1 

It  can  be  obtained  that 

A"  =  i  D.fi  i  =  1,2.  ...6  (21) 

where  K  =  D . f . ,  D.  =  f.A..,and  f.  =  — 
ii  i  jii  l  . 

J  i 

Substituting  Eq.  (16),  (19)  and  (21)  into  Eq.  (15),  one  obtains  the 
constitutive  equation 


o  _  r  o 
o.  =  C ,  .  e. 

i  ij  3 


j  1 . 2  .... ,6 


(22) 


where  C^.  is  the  elastic-plastic  matrix  which  should  be  updated  for 
nonlinear  analysis  at  each  step. 

The  nonlinear  analysis  is  carried  out  by  the  incremental  method 
combined  with  iterative  correction  for  the  load  vector  at  every  load 
step  (see  [27]).  The  following  iterative  process  is  adopted: 

a.  the  load  increment  is  applied  and  the  elastic  stress  and 
strain  increments  are  obtained. 

b.  The  yield  function  f(o)  is  calculated  for  each  element.  For 
the  elements  where  f  '  0,  the  elasto-plast ic  strain  increments 
are  calculated  using  the  elastic-plastic  matrix  (see  Eq.(22)]. 
The  stresses  are  modified  in  the  proportion  of  the  strain 


increment  to  artificially  make  f  >  0.  The  d  i  if  erences  in  the 
stresses  are  treated  as  the  initial  stresses, 
c .  The  residual  load  vector  is  computed  from  the  initial  stresses, 
and  step  (b)  is  repeated.  The  iteration  is  carried  out 
the  residual  vector  becomes  very  small  (say  one  percent  of  the 
Load  vector  in  the  preceding  iteration).  Once  convergence 
is  achieved,  the  next  load  increment  is  added  and  the  whole 
iteration  procedure  is  repeated. 

RESULTS 

In  order  to  validate  the  elastic-plastic  model  used  here,  the 

problem  of  an  isotropic  plate  with  a  notch  (see  Fig.  2)  is  solved  by 

using  the  formulation  developed  here  for  the  three-dimensional  analysis. 

The  notch  stress  vs.  displacement  is  shown  in  Fig.  1.  The  yie Id  strength 

2  5  2 

of  the  material  is  30  kg/mm"  and  E  =0.2x10  kg/mm  and  v  =  0.3.  The 
present  three-dimensional  analysis  consists  of  72  brick  elements;  the 
maximum  notch  stress  at  the  end  of  six,  seven,  and  eight  iterations  is, 

respectively,  19,  21  and  23  kg/mm".  The  convergence  trend  indicates 

2 

that  the  maximum  stress  value  will  he  lower  than  23  kg/mm'.  Results  of 
a  two-dimensional  (plane  stress)  analysis  of  the  same  problem  is 

2 

available  in  Reference  23.  The  maximum  notch  stress  was  19  kg/mm  and 
this  was  obtained  by  using  245  two-dimensional  plane  elements  and  51 

iterations.  From  the  classical  slip  line  solution  of  this  problem  the 

2 

maximum  notch  stress  is  18  kg/mm  .  Considering  the  fact  that  coarse 
mesh  and  fewer  iterations  (8  iterations)  are  used  in  the  present  study, 
and  the  present  study  is  based  on  the  three-dimensional  formulation  of 
the  problem,  one  can  conclude  that  the  results  obtained  here  are 


reasonably  good. 


1 1 


A  symmetric'  three-layer  square  laminate  (0°/90°/0°)  of  dimen¬ 
sions  a  x  a  x  h  and  subjected  to  uni  form]  v  distributed  vertical  load 
of  intensity  on  the  top  surface  is  analyzed.  The  plate  is  simply 
supported  on  all  the  edges.  Taking  advantage  of  the  symmetry,  only  a 
quarter  plate  is  analyzed. 

The  material  of  the  lamina  is  assumed  to  have  the  following 
properties  with  respect  to  the  material  symmetry  axes: 

I. inear  Elastic  Parameters 


model  with  the  yield  stresses  X  =  10,000  psi  Y  *  4000  psi  and  Z  *  4000 
psi  and  2)  low  strength  model  with  yield  stresses  X  *  5000  psi,  X  *  2000 
psi  and  X  -  2000  psi. 


Plots  of  the  notui  imens iona 1 i zed  deflection  versus  nondimens ion- 
alized  load  are  presented  in  Figs.  1-  S  tot  s ide- t n- i h i ckness  ratios 
a/h  -  r) ,  10,  20  and  100.  The  figures  also  contain  results  obtained  by 
the  use  of  the  modified  Ramberg-Osgood  model  in  the  two-dimensional 
laminated  plate  theory  (see  [b,  13]).  For  a/h  ratios  5  and  10  (see 
Figs,  i  and  4) *  the  results  are  in  close  agreement  in  the  initial  stages 
of  the  load,  and  the  deflections  predicted  hv  the  plate  theorv  are  lower 
than  those  predicted  bv  three-dimensional  theorv  in  higher  load  levels. 
For  (a/h)  =  20  (see  Fig.  5),  the  results  obtained  bv  the  plate  theory  do 
not  agree  closely  with  those  of  the  three-dimensional  elasticity  theory. 
This  is  due  to  the  fact  that  full  integration  (2x2x?)  is  used  in  the  3-D 
analysis,  which  results  in  so-called  locking  o4  element  stiffness 
coef  f icients  [13]. 

The  deflections  obtained  in  three-dimensional  analysis  using  the 

Ramberg-Osgood  model  and  elastic  plastic  model  are  also  compared  in  Fig. 

4  for  a/b  =  10.  A  similar  comparison  is  shown  for  (a/h)  =  100  in  Fig.  6. 

It  is  observed  that  the  Ramberg-Osgood  model  is  able  to  predict  an 

average  of  the  deflections  produced  bv  the  plasticity  models.  Figures  6 

and  7  contain  similar  results  for  stresses  (  >  at  the  (enter  ot  the 

x 

plate.  Here  also,  the  Ramberg-Osgood  model  predicts  an  average  of  the 
stresses  compared  to  the  stresses  predicted  by  the  elastic-plastic 
models.  This  observation  is  encouraging  to  recommend  the  use  of 
Ramberg-Osgood  model  for  design  purposes  because  the  Ramberg-Osgood 
model  is  much  simpler  than  the  elastic-plastic  analysis  which  requires 
many  iterations  at  each  load  level.  These  figures  clearly  demonstrate 
the  effect  of  material  nonlinearity  on  the  behavior  of  a  simply 
supported  composite  plate  subjected  to  uniformly  distributed  load. 
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Figure  5.  Load-Deflec  tion  Curve  for  Cross -PI v  (0°/90°/0 
Equal  Phiekness  Lavers)  Square*  Plate  Under 
Uniformly  Distributed  Load  q  (S  =  a/h  -  20). 


o  Ramberg -Osgood  (  M)) 


Load  o  =  x  10 
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Figure  6.  Load  Deflection  for  Cross-Ply  (0°/90°/0°>  Equal 
Thickness  Lavers)  Square  Plate  Under  llniformlv 
Distributed  Load  q  (S  =  a/h  =  100). 
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igure  7.  Normal  Stress  vs.  Load  for  Cross-Plv  (0°/90°/0' 
Equal  Thickness  Layers)  Square  Plate  Under 
Uniformly  Distributed  Load  q  (S  -  n/h  =  10). 
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CONCLUSIONS 

A  fully  three-dimensional  analysis  of  cross-ply  laminated  plates 
with  geometric  nonlinearity  is  presented.  Two  types  of  material  nonlinear 
models  are  used:  (1)  uncoupled  one-dimensional  modified  Ramberg-Osgood 
relation,  and  (2)  elastic-plastic  model  with  modified  Hill’s  criterion 
for  anisotropic  media.  It  is  observed  that  the  uncoupled  Ramberg-Osgood 
model  was  able  to  predict  the  average  material  behavior  compared  to  the 
elastic-plastic  model.  The  df lections  obtained  by  the  three-dimensional 
analysis  using  the  uncoupled  Ramberg-Osgood  model  are  compared  with  those 
obtained  by  the  shear  deformation  plate  theory.  Reasonable  agreements 
are  seen  at  lower  load  levels  up  to  the  plate  side  to  thickness  ratios 
a/h  =  10. 
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